Renormalization flow towards gravitational catalysis in the 3d Gross-Neveu model 
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Catalyzed symmetry breaking arises from a parametric enhancement of critical fluctuations in- 
dependently of the coupling strength. Symmetry-breaking fermionic long-range fluctuations exhibit 
such an enhancement on negatively curved spaces, as is known from mean-field studies. We study 
gravitational catalysis from the viewpoint of the functional renormalization group using the 3d 
Gross-Neveu model as a specific example. We observe gravitational catalysis towards a phase of 
broken discrete chiral symmetry both on a maximally symmetric (AdS) and on a purely spatially 
curved manifold for constant negative curvature (Lobachevsky plane). The resulting picture for 
gravitational catalysis obtained from the renormalization flow is closely related to that of magnetic 
catalysis. As an application, we estimate the curvature required for subcritical systems of finite 
length to acquire a gravitionally catalyzed gap. 



I. INTRODUCTION 

Mass gap generation and (chiral) symmetry breaking 
in relativistic fermionic systems can arise from a variety 
of mechanisms which are often related to certain cou- 
plings or interaction channels becoming dominant. An 
apparent counter-example is catalyzed symmetry break- 
ing, first studied in the context of magnetic catalysis [l|- 
Hl, where mass gap generation is triggered by the pres- 
ence of a magnetic field even for arbitrarily small values 
of the interaction strength. This phenomenon can be 
understood in various ways, the essence being that the 
long-range fluctuations driving the symmetry-breaking 
transitions are parametrically enhanced, see [J] for a re- 
cent review. More concretely, a magnetic field induces a 
fermionic fluctuation spectrum with Landau-levels con- 
taining a zero mode that leads to an enhancement of the 
density of states in the IR and to an effective dimen- 
sional reduction favoring symmetry breaking. Magnetic 
catalysis has found a rich variety of applications both in 
particle physics (chiral phases of QCD) 0-fLTl and con- 
densed matter physics [l3 - ll8| . 

Another simple picture for magnetic catalysis has re- 
cently been developed within the framework of the func- 
tional renormalization group in the context of the 3d 
Gross-Neveu model In line with the fact that 

symmetry-breaking phase transitions are often related 
to fixed points of renormalization group transformations, 
also magnetic catalysis can be related to the behavior of 
RG fixed points as a function of the magnetic field. This 
RG picture has already successfully been applied in the 
context of QCD [H. 

In the present work, we verify the underlying RG pic- 
ture of catalyzed symmetry breaking in the context of 
curved spacetimes. The fact that symmetry breaking and 
mass generation in fermionic systems can be influenced 
by ne gati ve curvature of the spacetime has been realized 
early |2l| , and is meanwhile reviewed in textbooks [22| . 
The phenomenon is typically studied at mean-field level 
and occurs in many different fermionic models [23T - l35j . 
In [HI, the similarity to magnetic catalysis was realized 
in terms of an effective dimensional reduction mechanism 



of the spectral properties of the Dirac operator. This jus- 
tifies the use of the terminology "gravitational catalysis" 

In fact, we find in the present work that the effective 
dimensional reduction and the corresponding enhance- 
ment of the density of states in the IR is directly related 
to the fixed point structure as identified below. From 
this RG viewpoint, symmetry breaking arises as a con- 
sequence of the fact that the coupling value required for 
criticality becomes arbitrarily small as a function of the 
curvature (the catalyzer). Hence, any finite value of the 
fermionic interactions ultimately becomes supercritical, 
typically driving the system towards the ordered phase. 

We investigate this RG mechanism within the simple 
3d Gross-Neveu model in the present work. We con- 
centrate on two different curved backgrounds with con- 
stant negative curvature: a maximally symmetric space- 
time (Anti de Sitter) and a purely spatially curved case 
(Lobachevsky plane). For both cases, mean-field stud- 
ies are already available, see [I?} and [H, [Hj], respec- 
tively. Whereas the former allows for an analytic treat- 
ment in terms of simple functions, the latter is poten- 
tially relevant for curved layered condensed matter sys- 
tems. For instance, the exitonic or anti-ferromagnetic in- 
stabilities in graphite and graphene have been associated 
with quantum phase transitions falling into the 3d Gross- 
Neveu universality class [HI, (39| . As catalyzed symmetry 
breaking is manifestly driven by the long-range modes, 
the RG analysis allows us to estimate the required cur- 
vature in relation to the length scale of the sample. 

This paper is organized as follows: in Sect. UH we 
briefly introduce the model and its formulation in curved 
space. Section HLT1 is devoted to an evaluation of the RG 
flow in its simplest formulation. The manifestation of 
gravitational catalysis is discussed in Sect. IIVI We es- 
timate the influence of finite probe length on gravita- 
tional catalysis in Sect. E]by means of a pseudo-critical 
coupling. Conclusions are drawn in Sect. IVII Relevant 
technical details arc deferred to the appendices. 



2 



II. GROSS-NEVEU MODEL IN CURVED SPACE 

We aim at investigating the 3d Gross-Neveu model 
in curved spacetime with metric g^ u (Greek indices run- 
ning from to 2) and signature (—,+,+) using func- 
tional RG methods. The microscopic action functional S 
at some ultraviolet (UV) scale A depends on the bare 
coupling constant Aa, the N Grassmann- valued fields 
tp = (tp 1 ) and the N conjugated fields ip = (V> 4 ), 



N 



N 



2N 



V^ + ^W) 2 



(1) 



where J = fd d Xy?—g is a shorthand for the integral over 

the d-dimensional spacetime and g = det g^ is the de- 
terminant of the spacetime metric. The differential oper- 
ator y = 7 M V M is composed from a set of d 7 -dimensional 
gamma matrices satisfying the Clifford algebra 



(2) 



where I is the identity in spinor space. In the present 
work, we use the irreducible representation, such that 
e? 7 = 2 specifics the dimension of the 7 matrices as well 
as the number of Dirac components of the fermions. For 
our computations with fermions in curved space, we have 
actually used the Weldon formalism [ijj which can be 
viewed as a generalization of the conventional vierbein 
formalism. However, a mapping to standard vierbein for- 
mulas is straightforward. The covariant derivative in the 
Weldon formalism reads 



(3) 



accounting for covariance with respect to the spinor and 
spacetime structure. Here, V ^ is the afhne spin connec- 
tion, implicitly defined by 



(i) o = v^ = V + r;; p 7 p + [r^], 

(ii) o = trr M , 



(4) 



with Y" the Christoffel symbolsQ In the present work, 
we are interested in a discrete "chiral" Z2 symmetry, 
where the nontrivial transformation is defined by [42j 



ip(x) — > —ip(—x), ip(x) — > tp(—x). 



(5) 



This symmetry acts simultaneously on all flavors. It 
can spontaneously be broken by a chiral condensate 



1 More generally, the spin connection T M can additionally accom- 
modate a U(l) gauge field A-fi |4 jl . In this case, the right-hand 
side of (ii) in J4]l would be given by trT^ = —id^qA^. where q 
is the charge under the (7(1) gauge group. In the present work, 
we ignore such a (7(1) gauge sector, setting = 0. 



(ipip) 7^ 0, which for finite interactions goes along with 
a mass gap generation. Incidentally, the 3d Gross-Neveu 
model actually has a much larger continuous U(iV) flavor 
symmetry also allowing for more complicated breaking 
patterns [H]0 



III. FERMIONIC RG FLOWS IN CURVED 
SPACE 

In the following, we use the functional renormaliza- 
tion group to compute the RG flow of the Gross-Neveu 
coupling function of the (negative) curvature. We 
employ the Wetterich equation [46|, describing the flow 
of a scale-dependent effective action functional Tk as a 
function of an IR regulator scale fc, 



d k r k [r,^} = -STr 



(It 



d k Rk 



(0) 



The effective average action I^ is related to the stan- 
dard generating functional for IPI correlation functions 
in the limit Y = Y k ^. Towards the UV cutoff k -> A, 
Y k approaches the microscopic action. The rcgulariza- 
tion is encoded in the regulator function R k , see below. 
For reviews of the functional RG adapted to the present 
context, see Refs. (47145^ . 

In the present work, we evaluate the flow within a 
rather simple approximation for the effective action. For 
this, we truncate the effective action to 



r k W,^} 



^ + 7^ WO 2 



(7) 



where the only scale dependence lies in the four fermion 
coupling X k . Furthermore, X k parametrically depends on 
the curvature of the background manifold. The IR regu- 
larization is ensured by a chirally symmetric regulator of 
the form 



(8) 
(9) 




2 In many 3d condensed matter systems where the Gross-Neveu 
model is considered as an effective theory, the low-energy degrees 
of freedom can be arranged into Nf 4-component Dirac spinors, 
corresponding to a reducible representation of the Dirac algebra. 
This reducible representation can be constructed from a suitable 
combination of 2-component spinors such that N = 2N{ in terms 
of the counting of fermions of the present work, see, e.g., |44| for 
a review. Note that the Gross-Neveu interaction term considered 
in this work ~ (i/>V) 2 corresponds to ~ (tp"/4,5ip) 2 in the reducible 
4-component notation of |4:| (or to ~ (i/'735i/ 1 ) 2 in the notation of 
even N. The critical properties of the discrete chiral 
transition are, however, identical to a 4-component Gross-Neveu 
model with a (ipip) 2 interaction as considered in l4 r l . 
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where the superscript T denotes transposition in Dirac 
space, and S(x, y) represents a spin-valued delta distri- 
bution, keeping track of the spinor or conjugatc-spinor 
transformation properties associated with the spacetime 
arguments. In App.fAl we briefly summarize our conven- 
tions. 

For practical computations, we use a Callan-Symanzik 
type regulator, that facilitates the use of proper time 
representations, 



r(x) 



1 



(10) 



Within our investigations we restrict ourselves to nega- 
tive curvature, giving rise to gravitational catalysis. It is 
intuitively clear, that positive curvature (e.g., a sphere) 
generically suppresses IR modes and thus reduces the 
density of states of low lying modes [U [22|. We also 
consider the system mainly in the large-TV limit. At fi- 
nite N , further pointlike fermionic self-interactions are 
generated which correspond to operators with an explicit 
curvature dependence^ We expect no qualitative modi- 
fications from these operators and hence ignore them in 
the following. This approximation becomes exact in the 
limit N — > oo. 

It is straightforward (cf. App. [Aj to calculate the flow 
of the coupling as an implicit functional of the choice of 
the manifold, which enters via the spectrum of the Dirac 
operator 



. 2\l 



STr 



(I + r)- 



S(x,y) 



(11) 



For this calculation, it suffices to project the flow onto 
constant fields ip l (x) = with d^ 1 = 0. Here, 

Q = [ 1 denotes the spacetime volume. The operator 

x 

occurring in the trace is related to the square of the reg- 
ularized fermionic Green's function in curved spacetime. 
This has a direct correspondence to a Feynman diagram 
representation, see Fig. [TJ as the flow in the present 
simple truncation is driven by a single fermion bubble 
(and RG-improved resummations thereof). In the fol- 
lowing we distinguish between the cases of a maximally 
symmetric spacetime in Sect. MI Al which can be treated 
fully analytically, and a negatively curved space in Sect. 
IIIIB1 which is a more interesting case in view of two- 
dimensional condensed matter systems. 




Figure 1. The diagram schematically exemplifies the fluctua- 
tion contributions to the flow of the Gross-Neveu coupling. 
The double lines represent the fermion propagator on the 
curved manifold. The full circles denote the RG-improved 
couplings, and the crossed circle marks the regulator inser- 
tion. 



side of Eq. (|11[) in terms of a simple proper time rep- 
resentation. Other shape functions would still permit to 
use a proper time representation but would lead to more 
intricate k dependencies. The Laplace transform of the 
operator of (TTTj) reads 



{1 + T) -2S(0 = _ k 4 



ds se 



(12) 



which upon insertion into (TTTJ) yields 



d k X k = i — 



ds s e 



STr 



s y2 5(x,y) 



-9 



(13) 



The expression inside the super trace is known as 
the heat kernel K(x,y;s) — c ls ^ S(x, y) / — g of the 
(squared) Dirac operator. It satisfies 



(i) d s K = iTK, 

S(x,y) 



(ii) 



lim K 

s\,0 



(14) 



-9 



and was calculated in [55} for any maximally symmetric 
space, with Euclidean signature. For our case the solu- 
tion to this equation can be obtained in an easier way 
with a special ansatz (cf. App. |B|) : 



A. Maximally symmetric spacetime 

The use of a Callan-Symanzik regulator shape func- 
tion, cf. Eq. (fT0|) . facilitates to rewrite the right-hand 



K 



cosh-i 



w 



\R\ 



s sinhio 12 cosh? 



(47^ 



-U, (15) 



where w 2 



\R\d, 



and (1g(x, y) is the geodesic distance 



between the points x and y. The parallel transporter U 
(Wegner- Wilson line) is defined by 



3 A similar mechanism has been observed in |19(I for the case of a 
magnetic field, and the corresponding operators have been clas- 
sified. 



U(x, y) = P exp 



, d^(i)„ 



it)) 



(16) 



4 



where P denotes the path ordering prescription and z(t) 
is the geodesic between x and y with z(t = 0) = x and 
z(t = l)=y. 
Now we are able to calculate the supertrace as 



This implies that the operator ^ 2 is separable into 



i S(x,y) 



STr 



and finally get 



2QN 



(4tt) 



s \s 12 



1^1 
24/c 2 



e 



(17) 



(18) 



(25) 



where the curvature is only induced by the spatial com- 
ponents. 

For a simpler calculation of the beta function, we per- 
form a Wick rotation x° — > —ix° and perform again a 
Laplace transformation of the operator occurring in Eq. 

(ED, 



In terms of the dimcnsionlcss coupling, 

Afc = k\ k , 
we obtain the beta function j3\ 



13 x = kd k \k = A fe 



A| 
2tt 



1 



\R\ \ 

24k 2 J 



(19) 



(20) 



This is an ordinary differential equation that paramet- 
rically depends on the curvature. It can be solved by 
straightforward integration: 



A A 



l-^(l-x) 



1 



\R\ 

24fcA 



(21) 



The initial value is given by the dimensionless coupling 
Aa which in terms of the initial bare Gross-Neveu cou- 
pling reads A a = AAa- 



B. Negatively curved space 

For the case of a manifold where the spatial part has 
a constant negative curvature, we choose a special set of 
coordinates such that the metric can be expressed as 




(22) 



where gij (Latin indices running from 1 to 2) represents 
the metric of a two dimensional maximally symmetric 
space and therefore only depends on the spatial coordi- 
nates. Hence the gamma matrices are time independent 
as well. The Christoffel symbols vanish for every time 
component = T ^ v = and also the curvature tensor 
vanishes if any index is zero R^ vp q = 0. From Eq. (|4]) 
(i), we infer that the spin connection Tq has to be pro- 
portional to I, 



= V07" = [To, 7" 



(23) 



Moreover, From Eq. ((4]) (ii), we conclude that Tq even 
has to vanish completely, 



trr = 0. 



(24) 



(I + r)- 



5(x,y) 



ds sc~ 



5(x,y) 



'9 J V~9 

o 

We arrive at the Euclidean analogue of Eq. (|T3 

,2 6(x,y) 



d k X k = ^ / 



-9 



(26) 



(27) 



Here, we calculate the super trace again with the aid of 
the heat kernel using that the differential operator ^ 2 is 
separable and the delta distribution factorizes in a time 
like and a spatial part 

esr s_{0 = c< (28) 

V^g Vg 

where g = det . The quantity x denotes the spatial 
coordinates of x and should be treated as a set of co- 
ordinates and not as a vector. Both factors satisfy a 
heat-kernel equation and can be solved analytically, cf. 



V°) = 



'4ns 



8(x° 

S(x,y) 
Vg " {4ns)iVW\ 



(29) 



2 cosh _1 f 



ve 2 =i«i cosh^ , 
dv =2=U, 
V cosh v — cosh w 



(30) 

where w 2 = ■ R ^ Q > U is the parallel transporter for the 
spatial part and g?g is the nonnegative spatial geodesic 
distance between x and y with (x,y) — <1q (x, y) — [x — 
y ) 2 - Plugging these relations into Eq. (f2"T|) gives 



k^k 



A| 
2tt 



3(a k ), a k 



3(a) = -jds 



, av 
-wcoth — . 
2 



(31) 
(32) 



The s integral is an integral representation of the modi- 
fied Bessel function of the second kind Kq(v) (56|. There- 
fore, we have 



3(a) = — dv vK (v) coth — . 

IT J 2 





(33) 
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With these results we are again able to derive the beta 
function for the dimensionless coupling Afc = k\ k , 

(3x = kd k Xk = \k-^-3(a k ). (34) 

The integration of this ordinary differential equation de- 
pending paramctrically on the curvature can be cast into 
an integral representation, 

k X A f\R\ 

OA 

This result is qualitatively similar to the maximally sym- 
metric case of Eq. (1211) . 

IV. GRAVITATIONAL CATALYSIS 

Let us now analyze the consequences of the RG flows 
for the long-range properties of the Gross-Neveu model. 
For both background manifolds, f3\ considered as a func- 
tion of Afc is a parabola where the prefactor of the 
quadratic part is scale and curvature dependent, see 
Fig. [21 For vanishing curvature, the /3\ function van- 
ishes at the two fixed points X k = (Gaufiian) and 
A»(i? = 0) = A cr = 2ir which corresponds to the well- 
known critical coupling of the Gross-Neveu model in flat 
space in this regularization scheme [45l . |54| . This criti- 
cal coupling separates the symmetric phase for Aa < A cr 
where the long range behavior is controlled by the non- 
interacting Gaufiian fixed point from the chiral symmetry 
broken phase for Aa > A cr - In the latter case, A& runs to 
large values towards the infrared. In the present simple 
truncation, X k in fact diverges at a finite scale &sb s i§~ 
naling the transition into the ordered regime. The scale 
kgB is thus characteristic for the physical scales in the 
ordered phase. In [57j it has been shown that /csb actu- 
ally agrees with the value of the dynamically generated 
fcrmion mass mt as obtained in mean-field approxima- 
tion. Since we are working in the large- TV limit anyway, 
we will use this mean-field identification in the following: 
TO f = k SB . 

The existence of the non-Gaufiian fixed point A cr can 
be attributed to the competition between the power- 
counting scaling (the linear coupling term in f3\) and the 
interaction terms ~ \\. In our RG picture, the interac- 
tion terms are enhanced by negative curvature as soon as 
the wavelength of the fluctuations becomes of the order 
of the curvature scale. As a consequence, the interacting 
second zero of the (3\ function no longer is a true fixed 
point but becomes scale dependent. This "pseudo-critical 
coupling" A p = A*(|i?|/fc 2 ) moves towards the Gaufiian 
fixed point for decreasing scale k, see Fig. [2] 

Any finite initial coupling strength Aa will eventu- 
ally become larger than A*(|-R| /k 2 ) for small RG scales 
k — > 0. By this mechanism, the system is forced into 
the symmetry-broken phase even at the weakest initial 



2.0 




1 2 3 4 5 



Figure 2. Plot of the RG f3\ function of the coupling \ k for 
different values of the scale dependent negative curvature 
(from top to bottom: 0; 5; 20; 100; 1000). The black lines de- 
pict the f3\ functions for the maximally symmetric spacetime 
(AdS), cf. Eq. (|20p : arrows indicate the flow towards the IR. 
The green dotted graphs show the flows for the case of purely 
spatial curvature (Lobachevsky plane), cf. Eq. p4[) . In addi- 
tion to the Gaufiian fixed point, there exists a non Gaufiian 
fixed point (full red circle at A cr = 2tt for the present regu- 
lator scheme), separating the symmetric phase for Aa < A cr 
from the broken phase for Aa > A cr for zero curvature. For 
finite curvature, this critical point becomes scale-dependent 
and moves towards the Gaufiian fixed point for increasing 
scale dependent curvature, i.e., with decreasing IR scale for 
fixed curvature. In the case of vanishing curvature, the sym- 
metry is preserved and no mass is generated for initial values 
Aa in the blue dashed region. 



coupling. We observe this mechanism in both cases of 
negative curvature, the maximally symmetric as well as 
the purely spatial curvature case. As we see in Fig. [5]thc 
influence of the curvature is somewhat stronger in the 
maximally symmetric case. 

As discussed above, we calculate the symmetry break- 
ing scale fcsB by searching for a zero of the inverse cou- 
pling. The fermion mass mf corresponding to this scale 
where the RG flow enters the symmetry-broken regime 
can thus be computed from the criterion 

A^ mr (|i?|/m f 2 ) = 0. (36) 

Upon partial bosonization (Hubbard-Stratonovich trans- 
formation), A^" 1 is related to the mass parameter of a 
composite bosonic field. Hence the divergence of the 
fermionic self-interaction simply corresponds to onset of 
the order parameter (54l . l58l , |59| | . Let us now analyze the 
two different backgrounds under consideration in detail. 



A. Maximally symmetric spacetime 

In the maximally symmetric case, the fermion mass 
defined by the criterion Eq. (|36p can be straightforwardly 
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computed from the running coupling Eq. (|21[) . yielding 



m,f 
T 



\R\ 
48A 2 




1^1 

48A 2 



1^1 

24A 2 ' 
(37) 



Plots of this gravitationally catalyzed fermion mass are 
shown as a function of the curvature as solid lines in 
Fig. [3] Let us discuss this result in various limits. In the 
zero-curvature limit R = 0, wc find m,f = for Aa < 2n. 
For super-critical couplings Aa > 27r, we rediscover the 
standard mean-field result in 3d, 



2n 

m f .o = rnf(i? = 0) = A ( 1 — — 

Aa 



(38) 



This is in perfect agreement with the known behavior in 
flat spacetime. 

Provided the fermion system is initially weakly cou- 
pled, Aa <C A cr = 2-7T, a leading order expansion can be 
performed for any value of the curvature, resulting in 



A 



mf 



48ttA 
\R\Xa 



for A A A < 1, 



(39) 



where we have reinserted the dimensionful initial cou- 
pling Aa = Aa/A. If we additionally consider the limit 
of small curvature, we find a linear dependence of the 
fermion mass on both the curvature as well as the cou- 
pling, 



the limit of large 



(40) 



By contrast, in the limit of large curvature, 
|i?|/(487rA 2 ) > 7r/A A > 1, wc find that m f A. 
In other words, large curvature induces immediate chiral 
symmetry breaking, such that the induced mass becomes 
of the order of the cutoff. Incidentally, this result is 
similar for the large-coupling limit: for Aa 3> 2tt, we 
again find that mf — > A to leading order independently 
of the curvature. 

The above results display explicit UV cutoff and 
regularization-scheme dependencies. Since the 3d Gross- 
Neveu model is asymptotically safe and thus non- 
perturbatively renormalizable [l9(, we can remove the 
UV cutoff by keeping an IR obvscrvable fixed while send- 
ing A — > oo. This "line of constant physics" defines a 
renormalized trajectory. This can most conveniently be 
done in the super-critical regime where Aa > 2ir such 
that the fermion mass in flat-space m^o of Eq. (j3"5)) de- 
fines a natural IR renormalization point In this case, 



4 In the sub-critical regime, the model is quasi conformal and can 
be renormalized, e.g., by fixing the coupling Aj. at a suitable 
renormalization point fc = /jtoa specific value. 



the generated fermion mass in the limit A 
written as 



co can be 




24m 2 



(41) 



where the first line holds for arbitrary curvature, and the 
second line represents a weak-curvature expansion being 
in perfect agreement with (2tJ ■ 

We emphasize that the fermions acquire a mass mf > 
for any given Aa as long as the curvature is nonvanish- 
ing. While we expect the tendency to drive the fermion 
system towards the broken phase through gravitational 
catalysis to remain also beyond our truncation, fluctua- 
tions of bosonic composites entering beyond the large- N 
limit typically provide for an opposite tendency. Hence, 
the status of gravitational catalysis beyond mean-field 
remains an interesting question. Analogously, the effects 
of beyond-mean-field fluctuations on mag netic catalysis 
arc under active current investigation |20L l6C| . 



B. Negatively curved space 

In the case of pure spatial curvature, the criterion 
Eq. (|36p cannot be resolved analytically, but we can give 
an implicit equation for the induced fermion mass mf 



2tt 
Aa" 



3(a) 



da, 




(42) 



which can be solved numerically. Though the basic pic- 
ture does not differ much from the maximally symmetric 
case, there are some interesting differences. As can al- 
ready be inferred from the beta function in Fig. [5] (dot- 
ted lines), the curvature induced mass in the spatially 
curved case is smaller than in the maximally symmetric 
case, cf. Fig. [3] Several limits can be discussed in an 
analytic fashion, using the series representation of 3(a) 
developed in App. [C] Furthermore, we need the integral 
5(a) defined by 



(i) 



3(a) 



a- 



da, 



5(a) - ^ ) 



(43) 



where (ii) fixes the constant of integration. The explicit 
calculation is done in App. [Cj The defining equation for 
the fermion mass in terms of £ is 



J(Om,) = ^(aA) + 



2tt 
Aaoa 



(44) 



This representation provides us with some interesting in- 
sight. First, we can show that there exists a unique 
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Figure 3. Gravitationally catalyzed fermion masses (mean- 
field level) as a function of negative curvature in units of 
the UV cutoff. The solid black lines display the maximally 
symmetric case, whereas the purely spatially curved case is 
shown as green dotted lines. The sets of three different lines 
correspond to super-critical, critical, and sub-critical bare 
fermion couplings, Aa — 1.6A cr , A cr , 0.8A cr from top to bottom 
(A cr = 2tt). As long as the background manifold is negatively 
curved, \R\ > 0, a finite fermion mass is generated. 



solution to this equation with < m; < A for any 
given negative curvature \R\ > and Aa > 0. This 
can be seen in two steps. The uniqueness is because 
the function $(a) € (-co, oo) for a€(0,oo) is bijective. 
This property holds, because 3(a) is positive and there- 
fore $ is strictly monotonically increasing. Owing to 
3(a) /a 2 — > l/a 2 for small a — > and 3(a) /a 2 — > l/(ira) 
for large a — > oo, the function 5" is not bounded. Since 
by assumption 27r/(AACtA) > 0, it follows from Eq. ([4"4")l 
that $(a m{ ) > v?(o!a) has to hold, which implies that 
ot mt > aA because of the monotonic behavior. This 
demonstrates that < mi < A as claimed above. 

Let us first check the flat spacetime limit \R\ — > 0. In 
complete agreement with Eq. (|38|) . we find again that 



^ = lim [a A vJ(a mf ) - (Xa-5(qla)] = 1 - —r^-, 
Aa \r\-k> A 



(45) 



where we have used Eq. (|C9|) . 

In the weak coupling regime, we expect the fermion 
mass to be small compared to the curvature scale, cf. 
Eq. (|40|) . Hence, we need $(a) for large argument as 
provided by Eq. (|CT0| . 



S a> 1) ~ . 

7T 



(46) 



Using this asymptotic behavior for 5(« m[ )i we infer from 
Eq. (031) that 




71 I 27r • 
t h a.A${a A 

Ot-A \Aa 



«A = 



\R\_ 

2A 2 ' 

(47) 



If in addition, the curvature is small, i.e., a a < 1, we 
can use the corresponding expansion of Eq. (|CJ9|) and re- 
place ciaS^a) — > —1- Taking differences arising from 
the Dirac representation into account, the exponential 
inverse-coupling dependence is in perfect agreement with 
the results of [331 ] . 

Despite the overall similarities to the maximally sym- 
metric case, we observe that the weak-coupling and weak- 
curvature limit of the case with pure spatial curvature 
shows some distinct differences. In particular, there is 
an exponential non-analytic dependence of the fermion 
mass on the coupling as well as on the curvature. 

This difference is reminiscent to magnetic catalysis in 
d = 2 + 1 and d = 3 + 1 Q, where the fermion gap is 
analytic in d = 2 + 1, but shows an essential singularity in 
the coupling in d = 3 + 1. Also in the present case, such a 
singularity shows up similar to BCS gap formation, as a 
consequence of the effective dimensional reduction of the 
fermionic fluctuation spectrum to d — >• 1 + 1 [36j. Also 
in this respect, our functional RG p icture is in perfect 
agreement with previous studies [33l 1371]. 



V. PSEUDO-CRITICAL COUPLING AND 
PROBE SIZE 

At zero curvature, the 3d Gross-Neveu model exhibits 
a critical coupling strength corresponding to a quantum 
critical point above which chiral symmetry is broken at 
large length scales. This critical coupling manifests itself 
as a non-Gaufiian fixed point of the RG flow. In the 
present regularization scheme, we identified A cr = 2n/A, 
or A cr = 2-7T in dimensionless conventions. As illustrated 
in Fig. [21 the fixed point strictly speaking no longer exists 
at finite negative curvature. The nontrivial zero of the /3\ 
function becomes scale dependent and eventually merges 
with the Gaufiian fixed point in the deep IR for k — > 
0, such that only the chirally broken branch of the f3\ 
function remains. 

Let us call the nontrivial zero of /3a a pseudo-critical 
coupling A p . We find 



A P = 



1 



\R\ 

2ik 2 



(48) 



for the maximally symmetric case, cf. Eq. (|20[) . and 

A 




(49) 



for the spatially curved case, cf. Eq. (j34|) . Since A p = 
A p (/c) is a monotonically decreasing function of scale k, 
the coupling Afc can eventually exceed A p such that the 
system becomes critical and runs towards the ordered 
phase. In this sense, 



Ap(fc c ) 



(50) 
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can be viewed as a criticality condition [6l|, denning a 
scale k c , where the system becomes critical. For lower 
scale, the system is driven towards the symmetry broken 
regime which is ultimately entered at ksB < k c defined 
above. The value of k c depends on the curvature as well 
as the initial coupling Aa (the latter is considered as ini- 
tially subcritical here and in the following). 

The preceding discussion implicitly assumed that k can 
run over all scales down to k = 0, such that the criticality 
condition ([5U|) can eventually always be satisfied. How- 
ever, if, for instance, the system has a finite volume char- 
acterized by a finite length scale L, also the fluctuation 
momenta are restricted, typically leading to an IR cut- 
off kh = 7r/£0 One may think of a finite probe length, 
such as, e.g., the size of a layer of graphene. This finite 
probe length L can lead to a screening of the gravita- 
tionally catalyzed ordered regime if k^ is larger than the 
would-be critical scale k c . Hence, A p (fcz,) can be viewed 
as a lower bound for the coupling strength required for 
symmetry breaking in a real system of finite length. It 
thus generalizes the critical coupling at infinite volume 
and zero curvature to the situation of finite volume and 
finite curvature. 

It is instructive to study A p as a function of the length 
scale L measured in units of a typical curvature length 
scale which we define by r = l/y/\R\. For finite sys- 
tems, this gives an estimate for how strongly a probe has 
to be curved in order to exhibit gravitational catalysis. 
In turn, for a given curvature of the probe, A p provides 
an estimate for the initial coupling strength required for 
symmetry breaking in the finite system. For instance, for 
the maximally symmetric case, we read off from Eq. (|48|) 
that 

\ p (k L =n/L)= \ 2 . (51) 

1 + 2471-^ 

For the spatially curved case, we incidentally find the 
same result in the limit of small curvature, i.e., L/r <C 1, 
using the expansion (|C3|) which is accurate even for val- 
ues of L/r~ 0(1). By contrast, the large curvature limit 
for the spatially curved case is different, cf. Eq. (|C7[) : 

A p (fc L = tt/L) = V2ir 2 ^\ CI , T - < 1. (52) 

From Eq. (|5T1) . it is obvious that probe length to cur- 
vature ratios up to L/r ~ C(l) lead to pseudo-critical 
couplings A p which deviate from the zero-curvature crit- 
ical coupling of the Gross-Neveu model only below the 
1% level. Significant deviations only occur for L being 
an order of magnitude larger than the curvature scale 
r. From the viewpoint of curvature-deformed condensed 
matter systems, a ratio of L/r ~ 0(1) appears to be 
"large" in the sense that - loosely speaking - a spatially 



5 Here, we tacitly assume that the boundary conditions are such 
that zero modes do not occur. 



curved 2d planar probe embedded in 3d Euclidean space 
would rather look like a 3d object. 

Another way to interpret these results is the follow- 
ing: consider a finite-probe system with a subcritical bare 
coupling, Aa < A cr , thus being in the symmetric (e.g. 
semimctal) phase. In order to gravitationally catalyze a 
transition to a broken (gapped or insulating) phase, the 
criticality condition (|50| has to be met for a sufficiently 
large k c > k^. In view of Eq. (|51[) this requires compar- 
atively strong curvature, i.e., a small curvature length 
scale r compared with the probe length L. 



VI. CONCLUSION AND OUTLOOK 

We have investigated the phenomenon of gravitational 
catalysis in the 3d Gross-Neveu model on specific mani- 
folds with constant negative curvature. While the mech- 
anism had already been studied frequently with mean- 
field methods as well as from the viewpoint of the fluc- 
tuation spectrum of the Dirac operator, we have added a 
new renormalization group picture to the comprehensive 
understanding of this phenomenon. The essence of this 
picture is that the critical coupling of the fermionic sys- 
tem, corresponding to a quantum phase transition in flat 
spacetime, is transmuted into a scale-dependent pseudo- 
critical coupling that flows to zero as a consequence of 
long- wavelength fluctuations (compared to the curvature 
scale). In this manner, the infinite- volume fermion sys- 
tem becomes critical for any arbitrarily weak coupling. 

We have identified the RG (pseudo-) fixed point mech- 
anism for two example manifolds of constant negative 
curvature: a maximally symmetric spacetime (Anti de 
Sitter) and a purely spatially curved case (Lobachevsky 
plane) . Both manifolds support the mechanism of gravi- 
tational catalysis, but exhibit a rather different behavior 
as far as the dependence of chiral symmetry breaking 
on the coupling and the curvature are concerned. The 
maximally symmetric case shows a linear dependence (to 
leading order) on both quantities which makes clear that 
the phenomenon is essentially perturbative, being remi- 
niscent of the "quantum anomaly" for fermions in a mag- 
netic field [62]. By contrast, order parameters indicating 
the symmetry broken state such as the induced fermion 
mass exhibit an essential singularity in both the coupling 
and the curvature for the case of the purely spatially 
curved case. This is in many respects similar to BCS- 
type gap formation. Again, our renormalization group 
picture goes hand in hand here with properties of the 
fermionic fluctuation spectrum, as analyzed in (36j . 

As a benefit, the functional renormalization group also 
gives a simple access to systems of finite extent by identi- 
fying the RG infrared cutoff with an inverse length scale. 
In this manner, we can estimate the fate of gravitational 
catalysis in finite systems. In fact, the phenomenon only 
occurs, as long as the curvature radius is sufficiently small 
compared to the systems length scale. We have been able 
to phrase this statement quantitatively by introducing a 
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pseudo-critical coupling. Thinking in terms of curved lay- 
ered condensed matter systems, rather large curvatures 
are needed compared to a realizable probe length in or- 
der to drive a sub-critical system into a phase dominated 
by gravitational catalysis. 

We would like to emphasize that an immediate appli- 
cation of our results to condensed matter systems would 
only be possible for reparametrization invariant systems 
such as fluid membranes [63j , curvature effects of which 
can be mapped onto the language of Riemannian ge- 
ometry. For tethered membranes or general lattice sys- 
tems, further phenomena connected to extrinsic curva- 
ture or curvature related defects can be become relevant. 
In this context, it is interesting to note that an exter- 
nal strain exerted on a graphene sheet in flat space in- 
duces a pseudo- magnetic field [64| , that may also support 
(pseudo-)magnetic catalysis. For such systems, we hence 
expect an interesting interplay between these various ef- 
fects if we expose them to negative curvature inducing 
strain. 

Finally, gravitational catalysis may become relevant in 
the context of asymptotically safe quantum gravity [65j . 
In conjunction with fermionic degrees of freedom [66[, the 
UV fixed point determining the shape of the universe at 
highest energies might go along with a negative (though 
scale-dependent) curvature (67j |. Whether or not gravita- 
tional catalysis in connection with gravitionally modified 
critical fermion interactions [H[ could become active and 
impose constraints on the matter content of the universe 
then is a highly involved question that deserves to be 
investigated in greater depth. 
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This gives rise to the following properties of the building 
blocks of the Gross-Nevcu action: 



(i) V^ = d^-VT M , 

(ii) (^r = ^v, 



(ni) 



(A3) 



Choosing the spin metric h to be anti-hermitean in 
Eq. (|A2[) (i) is convenient for our metric and Clifford 
algebra conventions, but differs from [4l|. Also, the ab- 
sence of any imaginary factor of "i" in the fermion kinetic 
term is due to these conventions. 

The flow equation © uses a rather condensed nota- 
tion. More concretely, we work in field space parameter- 
ized by the collective fields 



(A4) 



representing Grassmann-valued functions on the man- 
ifold, reminiscent to Nambu-Gorkov spinors. For in- 
stance, the classical action amended by the regulator 
term which is used for derivin g th e flow equation from 
the functional integral (cf. |47H54| ). reads in these con- 
ventions 



S k [d ) } = S^} + l 



</>{x)Rk{x,y)(j)(y). 



(A5) 



These conventions differ slightly from those commonly 
used, see, e.g., (59j . but turn out to be advantageous for 
coordinate space computations on curved manifolds. For 
instance, the representation of the unit element in field 
space becomes rather intuitive, 
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t{x,y) = 0(x) 



8<f>(y) 8<j>(x) 



0(2/) = 




By 8(x,y), we denote the spin-valued delta distribution, 
which fulfills 



Appendix A: Derivation of the beta function 

Let us first introduce our conventions. In the Weldon 
formalism [ill ], the Dirac conjugate spinor %jj is related 
to the hermitean conjugate of the spinor tp via the spin 
metric h 



a) = 1 6 -^dm, 

Jy v~9 

(n) fa) = ( m 5 -^-- 

Jy V-9 



(A7) 



Therefore, 



which is implicitly defined by 



(Al) 



(i) h) 

(ii) it 
(iii) 



-h, 

-hj^h' 1 , 
■ dfth - hT^ 



(A2) 



V (*)fe/) = gi^jElg + T ^ S A?M (A8) 

has to hold. With this notation, we indicate that S(x, y) 
transforms as a spinor in x and a Dirac-conjugated 
spinor in y, i.e. S(x,y) can be interpreted as 8(x,y) = 
8{x — y)U{x 1 y), with the standard scalar delta distribu- 
tion 8 (x ~ y). We remark that there is a difference be- 
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tween 5(x,y) and S(y,x) T in the spinor structure, since 



Appendix B: Derivation of the heat kernel 



I" 



,S(y,x) T _ 9 ( x) S(y,x) T rT 6(y,x) T 



v (x)%>^) 



(A9) 



(A10) 



For the evaluation of the flow equation Eq. (O, we pro- 
ceed in a standard fashion. We decompose 



If ' + R k = T k + v k 



(All) 



into a field-dependent part T k and a field-independent 
part Vk in order to expand the Wetterich equation in 
powers of the fields, 



Following [55[ , we choose for the heat kernel satisfying 
Eq. (Ti"4")) the ansatz 

K(x, y; s) = f(d G (x, y), s) ■ U(x, y) (Bl) 

where / is a scalar function of d G the geodesic distance 
and the proper time s. Plugging this into equation (|14[) 

using A 



— coth(2iu), B = J ^ tanhu) and 



(i) = 1, = dftdc, 

(ii) d^riy - T^rip = A(g tlv - n M n„), (B2) 

(iii) V fi U = B[ lfl , lv ]n v U, 



5 fc r fc = i £(-!)" STr [(n -1 ^^)^" 1 ^)"] ■ ( Al2 ) 



n=0 



Since 



A*/-[(M)I + *tt] 



V k = ifc 



v^(l + r(r)) 











T (l + r(r T )) 



1, (A14) 



we observe that only the term n = 2 in Eq. (|A12[) can 
contribute to the flow of X k - We need 



/ rr (t) q \ 



(A15) 



with r'(x) = -£-r(;r), and therefore only the diagonal of 
(■p^T 1 ^) 2 is required. In the limit N — > oo, only the 
following terms remain: 



[(V^F k f] ll {x,y) = - 



-(l + r(r)) z V-9 
[{V^T k )\ 2 {x,y) = [{V^T k )%{y,x). 



The LHS of Eq. (|aT2)| boils down to 
d k T k 

and the RHS yields 



2iV V 7 ' 



(A16) 
(A17) 

(A18) 



iSTi:^- 1 ^)^- 1 ^) 2 ] 



2A 2 

i^-^f^^) 2 STr 

A 2 fc 3 ^ J 



r'(r) <5(x,y) 
_(l + r(r)) 3 7=3 



(A19) 



Inserting the Callan-Symanzik regulator (|10p , we end up 
with Eq. (TTTI) of the main text. 



cf. [55|, we get 



R\ 



Q=(^-8B 2 )f + 2Af + f" + if (B3) 



with /' = dd G f, f = d s f. Because of the boundary 
condition for K and the regularity of U(x, y = x) = I we 
observe that 



lim f 



5(x - y) 



-fJ 



(B4) 



has to hold, with S(x — y) representing the standard delta 
distribution. One suitable representation for the delta 
distribution in 3d maximally symmetric space with neg- 
ative curvature in the limit s \ is 



S(x — y) e 1 4 
— 7=^ = Ta — ex P I 1 a 

V~9 (47TS)2 



(B5) 



Next, we factorize / into the delta part and an auxiliary 
function p of d G and s, 



f{d G ,s) = 7=— 7— r rexp 1— 

V s (47r)5 V 4s 



(B6) 



Expanding p in powers of ^ and using the boundary con- 
ditions, 



leads to 
(i)0 



p = -pi(<h) + m(d G ), 

s 



h- A ) Pi+p[ ' 



(B7) 



(ii) = d G Ap + d G p' + [8B 2 -^) Pl -2Ap' 1 -p'l, 



(m)0^^8B 2 -^y o -2Ap / -pl 



(B8) 
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with the boundary condition pi(0) = 1. From equation 
(i) and the boundary condition, we see 



A similar approximation can be done for large a by 
expanding 



2w 



Pi 



sinh(2u;) 
Plugging this into (ii) gives 

\R\ 



(B9) 



Po 



12 cosh" w 



(BIO) 



where we have eliminated the constant of integration us- 
ing equation (hi). From this, we finally get 



K 



e' 4 3 



w . . \R\ 
l- 



coshw \s sinhui 12 cosh wj (Aft)zyfs 



e '4 



-U. 



(Bll) 



Appendix C: Curvature expansion 

For the detailed analysis of the spatially curved 
case, the integral representation of the running coupling 
Eq. (|35[) can be studied in various limits analytically. 
More specifically, 3(a) as defined in Eq. (j3"3")l can be ex- 
panded for small and for large values of a. Let us first 
consider an expansion of this function about a = 0, start- 
ing with an expansion of the integrand, 



) coth ■ 



2+l(avf-^(av) i + 0((av) 6 ). 



(CI) 



Using the standard integral 
dvv 2k K (v) = 2 2k - 1 T 





2 = 77 






22fc+l 





(C2) 

the small a expansion of 3(a) can be computed to any 
order. To order a 4 , we find 



9 4 

3 ^ 1+ 12-80 • 



(C3) 



Due to the factorial growth of the coefficients, cf. Eq. 
(|C2[) . the expansion is an asymptotic series. By compar- 
ison with the numerical result, the accuracy of Eq. (|C3j) 
turns out to be above 99% up to a ~ 1. 



coth 



a 1 ! 



1 



- = 1 + 2V ( 



2 1-e 

n— 1 

which holds for any av > 0. Next, we employ 

oo 

— 1 na arcosh(na) 



(C4) 



dvvK (v)e- no 



(™) 2 -l (( na )2_i)i 



(C5) 

which can be used for na > — 1. For —1 < na < 
1, an analytic continuation into the complex is im- 
plicitly understood, leading to a replacement of the 
term arcosh(na) ((na) 2 — l) 3//2 by [— arccos(na)] (l — 
\2\- 3 / 2 



(na) ) "'" here and in the following. This leads to the 
convergent series 



. a I sr^ 
3 a = - + - > 



(na) 2 



na 2 arcosh(na) 
((na) 2 — 



(C6) 



for any a > 0. Neglecting orders higher than i and using 
arcosh(na) — > ln(2na) for na — > oo, we arrive at 



_ , , a 7r In a n ( , A 

3(a) ^- + ^ + - l+7-ln — 

7r 3 a 3 V ir 



■ 1 

! 

Q 



(C7) 



where 7 w 0.577 is the Euler-Mascheroni constant and 
A w 1.282 is the Glaisher-Kinkelin constant. The accu- 
racy of this result is above 99% for a > 5. 

With the series (|C6|) it is possible to find a series rep- 
resentation for the required integral (j43| 



In a 

7T 



ln(2na) 



na arcosh(na) 



1 



(C8) 



This scries is convergent for any a > 0. The expansions 
for small and large a are 



a -> : ^(a) 
a — > 00 : T?(a) 



1 a a J 

"a +C+ 12 ~ 240' 
In a 7r In a 



7r 6 a 2 
where c is defined as the limit 

c = lim I #(a) + - ) w 0.364. 
a^o V a 



(C9) 
(C10) 

(Gil) 
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